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Following the identification of the tt bond in graphene, in this work, a tt bond constructed through 
side-to-side overlap of half-filled 6pz orbitals was observed in a non-carbon crystal of Pbi_a,Sna,Se 
(x~0.34) (PSS), a prototype topological crystalline insulator (TCI) and thermoelectric material 
with a high figure-of-merit {ZT). PSS compounds with a rock-salt type cubic crystal structure was 
found to consist of a bond connected covalent chains of Pb(Sn)-Se with an additional tt bond that is 
shared as a conjugated system among the four nearest neighbor Pb pairs in square symmetry within 
all {001} monoatomic layers per cubic unit cell. The tt bond formed with half-filled 6pz orbitals 
between Pb atoms is consistent with the calculated results from quantum chemistry. The presence of 
TT bonds was identified and verified with electron energy-loss spectroscopy (EELS) through plasmonic 
excitations and electron density (ED) mapping via an inverse Eourier transform of X-ray diffraction. 

PACS numbers: 61.50.Ah; 71.15.Ap; 79.20.Uv; 87.59.-e 


In search of high ZT compounds for thermoelectric 
applications, an intricate balance among the electric con¬ 
ductivity, thermal conductivity, and Seebeck coefficient 
is expected, which often leads to materials falling in a 
critical range between narrow band gap semiconductors 
and semimetals. [l[ The latest findings on topological in¬ 
sulator (TI) and topological crystalline insulator (TCI) 
materials show that these also fall into nearly the same 
critical range. [H Extensive studies of TI and TCI ma¬ 
terials have mostly investigated symmetry, and analy¬ 
sis has been performed from the perspective of energy- 
momentum in reciprocal space. Q A more intuitive in¬ 
terpretation from a real-space chemical bond model that 
shows the actual electron density distribution is desirable 
and would be especially useful for materials scientists and 
chemists focusing on materials design for future high ZT 
thermoelectric and TCI materials. 

PbTe-based materials exhibit a demonstrated high 
thermoelectric figure-of-merit {ZT) at high tempera¬ 
tures and are therefore good candidates for thermoelec¬ 
tric device applications. [l|, ^ Recently, Pbi-a^Sn^^Se and 
Pbi_a^Sna,Te with cubic symmetry have also been pre¬ 
dicted and confirmed to be topological crystalline insu¬ 
lators (TCIs) for X within a critical narrow band gap 
range, where band inversion and spin-orbital coupling 
occur for the crystalline symmetry-protected topolog¬ 
ical orders. [1-0 However, for PSS compounds with a 
rock-salt type face-centered cubic (FCC) Bravais point 
lattice (Fig. [U^a)), several chemical and physical fea¬ 
tures have yet to be characterized. First, PSS com¬ 
pounds of rock-salt type in FCC symmetry do not eas¬ 
ily cleave along the close-packed {III} planes with the 
largest inter-planar distance and, instead, always cleave 
along the {001} planes. 0, 0, [l0 This observation sug¬ 
gests that strong anisotropic chemical bonding must exist 


despite both Pb/Sn and Te/Se atoms exhibiting undis¬ 
torted octahedral coordination within a cubic rock-salt 
structure. Most surprisingly, surface mirror symmetry 
breaking along the (110) directions has been observed in 
recent high-resolution STM studies of a PSS compound 
in the critical TCI regime at low temperature, 0 which 
creates a puzzling conflict with the persistent 3D cubic 
symmetry of the bulk. [i,[nl 

Graphene has been shown to possess a conjugated tt 
bond system and a Dirac-cone-shaped surface band. [HI, 
in In this study, we find that the tt bond can also be 
identified in the confirmed TCI material Pbi-^n^^Se 
(x^O.34) with the presence of Dirac cones.pl, l7| Di¬ 
rect correlation between the tt band and Dirac cone 
via pseudo-spin has been previously established for 
graphene, 0 which is consistent with the conjugated 
nature of the tt bond at the molecular orbital level. 0 
To establish the intimate relationship between the con¬ 
jugated TT bond system and the presence of the Dirac 
cone, especially when the Dirac cone implies a gapless 
semiconducting behavior, we must first demonstrate the 
existence of the conjugated tt bond in PSS. The elec¬ 
tronic structure of the PSS compound was explored using 
X-ray electron density mapping and electron energy-loss 
spectroscopy (EELS) in this study. A molecular orbital 
chemical bond model is proposed to explain the observed 
experimental results, where PSS compounds are pro¬ 
posed to consist of covalent chains via cr-bond-connected 
Pb(Sn)-Se through hybridized s and p (sp) orbitals in the 
valence shell; moreover, conjugated tt bond condensation 
similar to that observed in graphene is proposed to exist 
between Pb along the (110) diagonal directions. 

Electron density (ED) mapping has been used to 
determine the actual valence electron distribution be¬ 
tween bonded atoms in condensed matter. The real 
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FIG. 1. (color online) (a) Crystal structure and (b) electron 
density contour mapping of Pbi-ccSuccSe (x~0.34) in 3D. (c) 
A thick cut of ED for the PSS-{001} layer between z=0.33 and 
0.83. The circle-shaped electron clouds (sashed circle in or¬ 
ange) in the diagonal direction between Pb(Sn) atoms are pro¬ 
posed to arise from the tt bond electrons, (d) The correspond¬ 
ing molecular orbital model for the PSS-{001} monoatomic 
layer in (001) projection; details are shown in Fig [3jc)-(d). 
The four bar-shape electron clouds shown in (c) correspond 
to the electrons in <j bond between Pb(Sn)-Se (dashed oval in 
red) and the lone pair electrons of Se (dashed oval in green) 
shown in (d); however, these two are not distinguishable for 
PSS in cubic symmetry. 


space electron contour can be determined experimen¬ 
tally with the aid of an inverse Fourier transform of 
the reciprocal /c-space information obtained via X-ray 
diffraction. The theoretical foundation of applying the 
Fourier method in crystallography has been thoroughly 
documented. [l5| The bonding nature and distribution of 
electrons within chemical bonds can be visualized us¬ 
ing this technique. [m, [13 Here, the ED of Pbi-a^Sn^^Se 
(x^O.34) has been extracted from the inverse Fourier 
transform of powder X-ray diffraction data using a step 
size of 0.1 at room temperature, as illustrated in 
Fig. mb). 

A thick slice (z=0.33-0.83) of the ED for PSS is shown 
in Fig. m^)- addition to the expected ED from elec¬ 
trons in CF bonds between Pb(Sn) and Se along the (100)- 
directions in cubic symmetry, it is puzzling to observe 
non-zero ED with a circular shape in the diagonal (110) 
directions between Pb(Sn) atoms. We believe that these 
experimentally observed electron clouds are contributed 
by unexpected bonding electrons that link the four near¬ 
est neighbor Pb(Sn) atoms in square symmetry per {001} 
layer. To explore the character of such an unconventional 
chemical bond between Pb(Sn), EELS studies were per- 
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FIG. 2 . (color online) EELS spectrum of Pbi-ccSn^Se 
(x~0.34). The inset shows the Ueff integrated from the 
energy-loss function, which indicates that the two peaks near 
7 and 15 eV correspond to n 7 r:(n 7 r+nc 7 ) = l:(l+2) electrons per 
atom participating in the e-beam energy absorption from tt 
and TT+cr plasmons. The 20-22 eV weak peaks are due to the 
04 , 5 -edge excitation of Pb. 

formed. 

EELS can provide useful information regarding the 
effective number of electrons participating in the col¬ 
lective oscillations within chemical bonds through en¬ 
ergy absorption of the tt or a plasmons. EELS spec¬ 
tra for graphite and Bi 2 Se 3 have been compared to pro¬ 
vide evidence of the volume plasmons corresponding to 
the collective excitations of electrons within tt and a 
bonds. 0 The assignment of tt and a plasmons has 
been satisfactorily verified in graphite on the experi¬ 
mental values of the effective number of electrons per 
atom participating in the collective plasma oscillation, 
ne//(co’), as 77 . 7 ^:( 77 . 7 ^+ncr)=l:(1+3) for the two peaks in 
low-loss region of the EELS spectrum near 7 and 27 eV, 
respectively. 00 This finding can be interpreted as 
there being 1 and 3 electrons per carbon participating 
in the tt and a bonding, respectively; these numbers are 
consistent with the molecular orbital model of three a 
bonds and one tt bond that is shared by three equiva¬ 
lent nearest neighbor carbon pairs in a honeycomb lattice 
construction, as illustrated in Fig. [3}b). The randomly 
distributed tt bond in space and time can be viewed as a 
conjugated system, similar to those observed in conduc¬ 
tive polymers. [23 

The EELS spectrum of Pbi_a,Sna,Se (x^O.34) is shown 
in Fig. [21 In parallel to the successful interpretation of 
EELS spectra for graphite and Bi 2 Se 3 ,[l 8 |, [l9| the low- 
loss peaks near 7 and 15 eV can be assigned as corre¬ 
sponding to the energy absorption by the tt and tt+ct 
plasmons, respectively. 7r-plasmon has also been revealed 
in the EELS spectrum of nanocrystal PbSe, 0 and the 
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two weak peaks at ^ 20-22 eV are identified as the Pb 
04 , 5 -edge excitation from Pb 5d electrons within an oc¬ 
tahedral coordination. The effective number of electrons 
participating in the collective plasma oscillations of tt 
and TT-j-cr plasmons is determined by the energy-loss func¬ 
tion analysis to be ^1:3 corresponding to the two peaks 
near 7 and 15 eV, respectively, as shown in the inset of 
Fig. H Following the similar assignment of and Ua- 
for graphite, [m the effective number of electrons partic¬ 
ipating in the two EELS low-loss peaks for PSS near 7 
and 15 eV, ^1:3 (inset of Eig. [ 2 ]), can be interpreted as 
77 . 7 ^:( 77 . 7 ^-hncr)=l:( 1 + 2 ), i.e., there being 1 and 2 electrons 
per atom participating in the tt and a bonding, respec¬ 
tively. Clearly, the proposed existence of tt bond in a 
non-carbon system PSS seems surprising, which requires 
a proper molecular orbital model to support. 

The crystal structure of PSS compounds has a rock- 
salt type cubic symmetry as shown in Eig. [IJa). All 
of the atoms in PSS compounds can be described us¬ 
ing an ECC Bravais point lattice with the selection of 
the (Pn/Sn)-(Te/Se) unit as a di-atomic basis, which is 
crystallographically identical to the typical NaCl ionic 
crystal. Because both Pb/Sn and Te/Se atoms have oc¬ 
tahedral coordination, conventionally the crystal struc¬ 
ture is drawn assuming that chemical bonds exist for all 
six (Pb/Sn)-(Te/Se) pairs within each PbSce or SePbe 
octahedron. However, the electronegativity differences 
between Pb(2.33)/Sn(1.96) and Te(2.10)/Se(2.55) are 
rather small to suggest the existence of a polar cova¬ 
lent bond between Pb/Sn and Te/Sn based on Pauling’s 
rule of ionicity. [HI The assumption of six chemical bonds 
per Pb/Sn within an octahedral coordination faces two 
contradictory facts; the first one is that the four valence 
shell electrons per Pb/Sn (Pb=[Xe]4/^^5(i^^65^6p^) are 
insufficient to form six a bonds with the six neighbor¬ 
ing Te/Se (Se=[Kr]4(i^^5s^5p^); and the second one is 
that for Te/Se in the chalcogen family below oxygen, 
at most two electrons are needed for stabilization into 
a filled inert gas configuration; i.e., a coordination of two 
is expected. It is apparent that the conventional crystal 
structure drawn with chemical bonds to link all (Pb/Sn)- 
(Te/Se) in octahedral coordination could be heuristic and 
requires further examination. 

As a prototype IV series semiconductor, the chemi¬ 
cal bond of Si crystal has been satisfactorily described 
through hybridized sp^ orbitals for all Si atoms in tetra¬ 
hedral coordination. However, the same bonding the¬ 
ory failed to apply directly to the IV-VI compounds of 
rock-salt type cubic structure, such as PbSe-PbTe-SnTe, 
mostly because of the partially covalent polarity due to 
the small electronegativity difference and the unpaired p 
orbital electrons. [23| In parallel to the well-known molec¬ 
ular orbital model for graphene as shown in Eig.[3fa)-(b), 
we propose a molecular orbital model that could maxi¬ 
mize the number of chemical bond for PSS to reach the 
condensed matter ground state following the valence shell 
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FIG. 3. (color online) (a) The electronic configuration and 
(b) molecular orbital model of graphene, where the a bonds 
are formed with half-filled sp‘^ hybrid orbitals per carbon in 
the honeycomb structure, and each tt bond is shared by the 
three nearest neighbor C-C pairs in both space and time, (c) 
The electronic configurations of Pb/Sn and Te/Se and (d) the 
corresponding molecular orbital model for a representative 
{001} monoatomic layer of the PSS compound, a bonds are 
formed with the half-filled hybridized sp orbitals per Pb/Sn 
and Te/Se along any equivalent axis in cubic symmetry (as¬ 
sumed x-axis here), tt bonds are proposed to form within both 
the equivalent xy- and xz-planes in cubic symmetry, as shown 
between the half-filled pz orbitals of Pb(Sn) atoms when xy- 
plane is chosen, through side-to-side orbital overlap along the 
(110) diagonal directions. 


electron pair repulsion (VSEPR) theory rigorously. [24 1 
Under the requirement of a coordination number of two 
for Te/Se in the chalcogen family, it is natural to assume 
sp hybridization from the s and p orbitals in valence 
shell near the Eermi level for both Pb/Sn and Te/Se, 
which leads to the formation of a covalent bond through 
the head-to-head overlap of half-filled sp orbitals as a a 
bond, as illustrated in Eig. [3fc)-(d). These covalent a 
bonds form chains along one of the three equivalent crys¬ 
tal axes in cubic symmetry following nucleation. How¬ 
ever, the remaining two pi {i = x^y^z) orbitals for Te/Se 
are both filled as inactive lone paired electrons, whereas 
the two remaining pi orbitals for Pb/Sn are both half- 
filled and subject to additional chemical bonding. One 
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FIG. 4. (color online) The calculated contours of the elec¬ 
tron density distribution for two atoms at equilibrium inter¬ 
atomic distance (lower) and 1.5 times farther apart (upper), 
(a) Two half-hlled 2pz orbitals between C-C atoms, (b) two 
half-filled 6pz orbitals between Pb-Pb atoms. The equilib¬ 
rium interatomic distances between C-C of 1.4 A is acquired 
from graphene, and the equilibrium distance between Pb-Pb 
is acquired from the PSS crystal structure of ^a~4.3A. For 
presentation purpose, the charge densities of (a) C-C and (b) 
Pb-Pb are multiplied by factors of 5 and 100, respectively. 


of the most probable chemical bonds in addition to a 
bond formation could be a tt bond formed between the 
half-filled 6pz orbitals of Pb atoms along the diagonal di¬ 
rection of x^-plane, as illustrated in Fig. [3fd), where a 
bond chain is named along the x-axis within the x^-plane 
and the other equivalent choice of tt bond formation via 
6py pairing within the xz-plane is not shown. 

Based on the tt bond formation via side-to-side over¬ 
lap of half-filled 2pz orbitals in graphene, we decide to 
explore the possibility of tt bond formation via side-to- 
side orbital overlap of half-filled 6pz (see Fig. [3l^c)-(d)) 
orbitals between Pb. Starting from the calculation of the 
probability function (|4^(x, 2 ;)p) of the corresponding 

niT bond orbital wave functions z)=N{ip + 

V'n'io), where ipnio is an npz orbital function of the atom 
of interest, a comparative study of the electron density 
distribution between 2pz and 6pz pairs as a function of 
inter-atomic distance is presented for C-C and Pb-Pb 
pairs in Fig. 01 The two atoms are located along the 


x-axis and the Slater rule is adopted for each atom to 
construct the npz wave function. Fig. HJ^a) shows the 
calculated contours in the ^ = 0 plane section of the 
electron density of two carbon atoms with C-C at vari¬ 
able distances. It is apparent that a tt bond is formed by 
the two half-filled 2pz orbitals when the inter-carbon dis¬ 
tance is close to the average C-C bond length (^1.4 A) for 
graphene. 0 Interestingly, for the PSS crystal with equi¬ 
librium Pb-Pb distance of ^a^4.3 A, Fig. HJ^b) shows 
that a TT-bond-like electron distribution can also form be¬ 
tween the two half-filled 6pz orbitals of Pb atoms with 
an inter-atomic distance of ^4.3 A. While the 27t bond 
orbital has no node along the 2 ;-direction and Gtt has mul¬ 
tiple nodes with nonzero electron density near each nu¬ 
cleus, the outer shell of Gtt does form continuous electron 
density between the 6pz pair resembling the typical 27r 
bond orbital. 

Considering the two unpaired electrons in each hy¬ 
bridized sp orbital of Pb(Sn) and Se, and the one un¬ 
paired electron within {i = p, z) orbitals of Pb(Sn) (see 
Fig. [3]^c)-(d) when a bond is named along the x-axis), 
each {001} monoatomic layer has two unpaired electrons 
in sp orbital per Pb(Sn) and Se participating in the a 
bonding, and one unpaired electron in pz or py (depend¬ 
ing on the choice of xy- or xz-plane, respectively) per 
Pb(Sn) participating in the tt bonding, which is in per¬ 
fect agreement with the interpretation of EELS spectrum 
of PSS (Eig.[2]) showing that ria ^ 2 for a plasmon and 
^ 1 for TT plasmon. The satisfactory agreement between 
EELS and the proposed molecular orbital model points 
to an obvious selection of two a bonds per atom that link 
Pb(Sn)-Se and one tt bond per Pb(Sn) to be shared by 
the four nearest neighbor Pb(Sn) pairs in square sym¬ 
metry for all {001} monoatomic layers in PSS crystal of 
cubic symmetry, as illustrated in Eig. [3}d). 

In addition to the evidence of tt bond existence from 
EELS, we find that the electron clouds along the (110) 
directions between the Pb(Sn) atoms in ED mapping (see 
Eig. djc)) can also be assigned to the high probability of 
TT bond electrons populated in between Pb(Sn) atoms. 
However, it is expected that the tt bond can only be 
formed one at a time among the four choices in space 
and time without destroying the original square symme¬ 
try. Based on the homogeneous ED contours at the mid¬ 
point between Pb atoms along the (110) diagonal direc¬ 
tion, the sharing must be evenly distributed at the four 
diagonal sites per {001} layer in space and time. While 
the electron is not dividable when fractional charge is im¬ 
plied, it is suggested that the tt bond has a conjugated 
nature similar to that observed in graphene. 0 In fact, 
if the TT bond is simply randomly distributed among the 
four sites per {001} layer permanently, i.e., the tt bond is 
formed without sharing as a conjugated system in space 
and time but randomly and permanently distributed in 
space only, no constructive interference from the diffrac¬ 
tion of TT bond electrons is expected. 
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In summary, a molecular orbital model is proposed for 
Pbi_a,Sna,Se, which is verified experimentally to show one 
unusual symmetry-protected tt bond sharing mechanism. 
The identification of conjugated tt bond system in PSS 
compounds, similar to that observed in graphene as a 
gapless semiconductor with a Dirac cone shape surface 
band, suggests that the topological crystalline insulator 
Pbi_a;Sna,Se has a great potential on device application 
via Dirac cone engineering, especially with the chemical 
approach on the conjugated tt bond system control. 
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